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It is shown by Monte Carlo method that the finite size 
scaling (FSS) holds in the two dimensional random-coupled 
Ising ferromagnet. It is also demonstrated that the form of 
universal FSS function constructed via novel FSS scheme de- 
pends on the strength of the random coupling for strongly 
disordered cases. Monte Carlo measurements of thermody- 
namic (infinite volume limit) data of the correlation length 
(£) up to £ ~ 200 along with measurements of the fourth or- 
der cumulant ratio (Binder's ratio) at criticality are reported 
and analyzed in view of two competing scenarios. It is demon- 
strated that the data are almost exclusively consistent with 
the scenario of weak universality. 

PACS numbers: 75.40.Mg, 75.10.Nr, 05.50.Jk 



I. INTRODUCTION 

Random disorder plays important role in many inter- 
esting phenomena of condensed matter physics. Some 
examples are Kondo problem and metal insulator tran- 
sition driven by random disorder. The two-dimensional 
(2D) randomly disordered Ising ferromagnet is the sim- 
plest nontrivial statistical model including the effect of 
another type of fluctuation in addition to usual thermal 
fluctuation. By the random disorder is meant either a 
random site dilution or random- valued positive coupling. 
The effect of the combined fluctuations of the thermal 
and quenched random disorder on the critical behavior 
of the system has been the main subject of the studies. 
The two dimensional random coupling (or random bond) 
Ising ferromagnet is defined by the Hamiltonian 



H — ^ JijSiSj 



J i3 > 0, 



±1, 



(1) 



<y> 



where the sum is over all the links of a lattice. To be more 
specific, for every link on a given lattice one assigns ran- 
domly to Jij one of two different ferromagnetic coupling 
constants - one that of the pure system J and the other 
perturbed coupling J'. In this way single realization of 
random distribution of coupling is made. Due to its ran- 
dom nature, obviously, the realization of the disorder is 
uncorrelated. We are interested in the thermodynamic 
behavior averaged over infinitely many different realiza- 
tions of the distribution of the coupling constants. 

McCoy and WuEl made the first successful attempt for 
such a model, and proved that the specific heat (C v ) is 
non-divergent in the 2D Ising system with one-directional 



and correlated random bond disorder. According to Mc- 
Coy and Wu, the divergence of C v is caused by the coher- 
ence of all the bonds acting together, so destroying the 
coherence by introducing one or two directional bond dis- 
orders makes C^finite. 

Many authorsfra, however, regarded the non-diverging 
behavior of C v in the McCoy- Wu model as a characteris- 
tic of the one-dimensional correlated disorder. Especially, 
Shalaev, Shankar, and Ludwig (SSLp, based on their 
observation that the continuum limit of the 2D random 
bond disordered Ising model is a certain type of Gross- 
Neveu model were able to obtain analytical prediction 
on the critical behavior of the disordered system. In this 
theory the disorder is regarded as a weak perturbation to 
non-interacting free fermion system (that is equivalent to 
the 2D pure Ising system) , so that the theoretical predic- 
tion is supposed to hold for weak random disorder. The 
main theoretical prediction can be summerized as: 



£ - t- v \\nt\ v , v = l, v = l/2 
C v ~ t~ a \ In | In a = 0. 

e-\ 17=1/4, 



X 



(2) 
(3) 
(4) 



where t is the reduced temperature and x 1S thermody- 
namic magnetic susceptibility. 

Note that in the context of SSL's theory the critical be- 
havior is modified by the disorder only logarithmically, 
with no difference in the values of the critical exponents 
from those of the corresponding pure system. In the pres- 
ence of extremely weak disorder the system must be al- 
most identical to the pure system, and the corresponding 
critical behavior must be almost indistinguishable from 
that of the pure system. This means that the modifica- 
tion of the critical behavior in this case must take place 
in extremely narrow scaling regime near criticality only, 
with the remaining scaling regime exhibiting the critical 
behavior of the pure system. The conventional belief is 
that Eq.(||)-Eq.([l|) remains to be valid for strongly disor- 
dered cases — only the range of scaling region where they 
hold becomes broader with the disorder. Hence the the- 
ory generally implies the presence of crossover from the 
critical behavior of the pure Ising system to the asymp- 
totic behavior Eq.(|)-Eq.(@). 

The expressions reflecting the presence of the crossover 
thus read 



£ ~ t- v [l + C\ hx(t)\f, u = l, u = l/2 
C v ~t- a ln\l + Cln(t)\+C, 

x ~e-\ »?=i/4 



(5) 
(6) 
{<) 
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The coefficients C and C as a function of the strength of 
disorder cannot be determined theoretically. Neverthe- 
less, obviously, C increases with the strength of disorder, 
and for a given value of C Eq.Q reduces to its asymp- 
totic form only when t is very small. 

These is a more general argument, originally due to 
Harrisu, and later elaborated by semi-rigorous approachO, 
on the effect of random quenched disorder to critical be- 
havior. The argument can be formulated in the context 
of the renormalization group theory, and it was obtained 
that weak, uncorrelated random disorder is irrelevant if 

Dv > 2, (8) 

where v is the critical exponent of the correlation length 
of the pure system. In our case D=2, v = 1, the random 
disorder has been regarded to be marginal leading to a 
logarithmic correction to power-law critical singularity, 
being consistent with the prediction of SSL. 

However, especially in 2D, a presence of marginal per- 
turbation in certain case such as in Kosterlitz-Thouless 
type phase transition leads to continuously varying crit- 
ical exponents. A rigorously known example of contin- 
uously varying critical exponents of 7 and is, but with 
the ratio 7/1/ remains unchanged from the value of the 
pure Ising ferromagnet, takes place in some subspace of 
the coupling constants of the 2D Ashkin- Teller model. 
In this subspace, the model is equivalent to a certain an- 
nealed ipandom coupling Ising ferromagnet. Also, many 
authorsErQ obtained qualitatively different expressions for 
the critical behavior of C v from that given by SSL. Ac- 
cording to these expressions, C v is non-diverging at crit- 
icality as is the case in the McCoy-rWu, model. 

It is now numerically establishectnlil that the value of 
r\ does not depend on the strength of disorder in both 
2D random-coupled and site diluted Ising ferromagnets. 
Another generic feature obtained from various numeri- 
cal studiesErtL2l is that other critical exponents such as 
7 and v increase with the strength of disorder at least 
effectively. In most cases this effective increase was in- 
terpreted as the multiplicative logarithmic correction. 
Whether these effectively varying critical exponents are 
due to crossover effect as claimed in Ref . ( p|,^2[-|T^) or are 
indication of genuine novel critical behavior induced by 
random disorder as claimed in Ref.( |l^,|l^) has been a 
controversial issue. Numerically it is very difficult to dis- 
tinguish a pure power-law criticality with small critical 
exponent from a logarithmic singularity, unless thermo- 
dynamic data in extremely deep scaling regime are avail- 
able. Previously claimed evidence for the predictions of 
SSL in most cases were nothing but the consistency with 
the logarithmic correction. This claimed consistency can- 
not rule out pure power-law singularity without the mul- 
tiplicative logarithmic correction. Especially it was ob- 
served that the thermodynamic values of the specific heat 
developa_a non-diverging peak at a temperature larger 
than /3 C E3 for very strongly disordered random site di- 
luted Ising model. This contrasts clearly rfflth Eq.(||), 
but agrees with other theoretical predictiontffl. 



In this paper we attempt to clarify the unresolved issue 
of the 2D random bond disordered Ising ferromagnet by 
numerically determining the functional form of the uni- 
versal ESS-, function defined in the context of novel FSS 
schemeHll3. As long as the FSS is valid, the form of 
the FSS function is unique for a given universality class. 
Accordingly, if the system is indeed in the same universal- 
ity class irrespective of the strength of disorder then the 
form of each finite size scaling function must be identical 
for different strength of random disorder. We will also 
demonstrate that the thermodynamic data of £ in suffi- 
ciently deep scaling regime tend to scale steeper than the 
asymptotic expression of SSL, Eq.(|2|). In the sections to 
follow we elaborate on the FSS method employed here, 
give a detailed description of our Monte Carlo simula- 
tion, report our results, and finally conclude with some 
discussions. 



II. GENERALIZED FINITE SIZE SCALING 

The fundamental assumption of FSS theoryB is that 
Az,(i), the value of some physical variable A on a finite 
lattice of linear size L, can be expressed as 

A L (t)=L"/"f A {s{L,t)), s(L,t)=L/£(t) (9) 

for a thermodynamic quantity A which has a power-law 
critical singularity A(t) ~ t~ p with the reduced temper- 
ature t = {fi c — 0) I '/3 C . Eq.(||) is valid for the values of 
L and £(t) that are sufficiently large; otherwise, there 
should be corrections to FSS. 

Notice that using the critical form for £, £(t) ~ t~ v , 
one can rewrite the scaling variable s(L,t) as 

s(L,t) = (A(t)/LP/ v )-^?, (10) 

so that Eq.(^|) may be rewritten as 

A L (t) = A{t)T A {s{L,t)\ (11) 

where the relation between the scaling functions j 1 a and 
J- A is given by 

Fa{8) = sr/»f A (s). (12) 

For A = £, Eq.([ll|) shows that £x,(t)/i is just a function 
of £(t)/L and vice versa, and this leads to the relationllSl 

A L (t) = A(t)Q A (x(L,t)), (13) 

where x(L,t) = £i,{t)/L is the ratio of the correlation 
length on a finite lattice to the linear size of the lattice, 
and Q A (x) is given by 

Q A {x) = FaU^Hx))- (1 4 ) 

Using the same observation, it is trivial to obtainE! an- 
other equivalent form, 
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A bL (t) = A L (t)G A (x(L,t)), 



(15) 



where & is a scaling factor and Oa(x) is another scaling 
function. 

It is evident that given ] a one can determine the other 
two scaling functions from Eqs.(12) and (|l4|). It is well- 
known that the standard scaling function Ja is universal, 
so that all the other scaling functions, J- a, Qa, and (Ll 
should be universal as well. It has also been arguedO 
that a certain asymptotic form of f A {s) can be expressed 
in terms of the critical exponent 5; by fitting this func- 
tional form one can extract an estimate for the critical 
exponent. 

It is worth stressing that use of the scaling function Q 
rather than J- or f A would be more convenient in many 
cases, particularly because one does not need the ther- 
modynamic correlation length to define the former. Note 
that there is no explicit t dependence of the scaling vari- 
ables, so that knowledge of the critical temperature and 
v is not required, and that x becomes independent of L 
at criticality. This L independent value of x at critical- 
ity, x c , which characterizes a universality class for a given 
geometry, forms the upper bound of x. In other words, 
the scaling function Q is defined only over < x < x c . A 
priori, the two limits of the scaling function Q are known 
for a continuous phase transition: lim^^o Q{%) — > 1 and 
lim x _,. Xc Q —>■ 0, because Al converges to its thermody- 
namic value in the former case while A(t) diverges in the 
latter case with Al{£) finite. Qa{x) turns out to be a 
monotonically decreasing function of x for A — \ an d £■ 

It is important to realize that the knowledge of the 
scaling function Q near x ~ plays as relevant role as 
that near x ~ x c to the extraction of necessary infor- 
mation of the critical behavior in (deep) scaling region. 
It can be easily seen by noting that x(L,t) for a fixed 
temperature arbitrarily close to criticality can be made 
arbitrarily close to zero by simply choosing a value of L 
sufficiently large. 

Eqs.([li"|) and (|l^) do not include any critical expo- 
nents, so that one might conjecture that their validity can 
be extended to non-power-law singularities. Although a 
general proof of this conjecture is missing, for the two- 
dimensional (2D) O(N) (N > 2) spin model which ex- 
hibits an exponential critical singularity, LuscherE3 ob- 
tained an explicit expression for the inverse correlation 
length (mass gap), 



&(*)~£oo(i) 



~4ttF(x) 



/x 2 x (1 + O (log L/L 2 )) 



1/2 



(16) 



Eq.(|l6[) is reduced to Eq.(|l5|) for the sufficiently large 
values of L as for C(log L/L 2 ) to be safely neglected. 
Note that 0(\ogL/L 2 ) is comparable to statistical errors 
of typical Monte ftpJp simulations at L=20 already. It 
is thus postulatedBlia that Eq.(jT|) is a generalized FSS 
form that is valid irrespective of the functional form of 
critical singularity, 



Man y n ovel applications of FSS were obtained thanks 
to Eq. (|l3|) . The key observation is that the scaling vari- 
able x has no explicit temperature dependence. First, it 
can be used to check the validity of FSS itself for a given 
physical system: one numerically calculates the scaling 
function Q^(x) at two arbitrary temperatures close to 
criticality. If the two scaling function thus calculated 
"collapse" unto a single function then Eq.(13) must be 
valid. The demonstration of the validity of FSS for the 
2D and 3D pure Ising models in this context can be found 
in Ref.( |l9|). The numerical determination of the 
were carried out for a variety of other statistical mod- 
els such as 2D Heisenberg and XY modelala. This check 
is much more unambiguous than thc-standard data col- 
lapse method of FSS on the modeled since it does not 
involve any adjustable parameters such as the critical 
temperature and the critical exponent of the correlation 
length v. Second, can be used to extrapolate infinite 
volume limit (thermodynamic) values of various physical 
variables as first demonstrated in Ref.( |lq,|2l|). 



III. THE SIMULATION 



We consider binary distribution of Jy. Namely the 
value of Jij at a link < ij > is randomly distributed 
between two positive values J and J' with probability 
p and I,— r p respectively. For p = 1/2, the system is 
sclf-dualcJ with the self-dual point given by 



tanh(J) = exp(-2 J') 



(17) 



A self-dual point equals the critical point of a system, 
provided the system has only one critical point. We 
fix J = 1 and p = 1/2 without any loss of generality, 
and consider three different values of J', i.e., J' = 0.9, 
0.25, and 0.1. Accordingly the self-dual points (critical 
points) are given by (3 C = 0.4642819 . . ., 0.80705185 . . , 
and 1.10389523 ... for J' = 0.9, 0.25, and 0.1 respectively. 

Our raw-data for each J' are obtained by choosing a 
realization of random J', then ruruaing Monte Carlo simu- 
lations in single cluster algorithmc3 with periodic bound- 
ary conditions; for each realization, measurements were 
taken over 10 000 configurations each of which was sepa- 
rated by 2-15 single cluster updatings according to auto- 
correlation times. The procedure is then repeated for dif- 
ferent realizations of J'. The average over all the different 
realizations converges as the numbers of the random re- 
alization increase; basically this mean value of a physical 
quantity is something physically interesting. To achieve 
the necessary precision for our FSS scheme the numbers 
of the different realizations we used are approximately 
20-40, 150-250, and 300-1000 for J' = 0.9, 0.25, and 0.1 
respectively; yet, in general, the fluctuation among differ- 
ent realizations of the random disorder is more significant 
than the statistical error for a given realization. This was 
particularly the case for J' = 0.1. Nevertheless, the av- 
erage over different realizations obviously converges to a 
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physically meaningful value. Our quoted error bars in our 
numerics are obtained by standard jack-knife method, 
taking into consideration the variation with different re- 
alization only. The smallest and largest values of L we 
used to extract our thermodynamic values are 20 and 
400 respectively. For J' = 0.05 it turns out that data 
of very large £ are not necessary, yielding compelling re- 
sults with the data for £ < 50 already. Thus the FSS 
extrapolation method is not used for this J 1 and 250-500 
different realizations of distribution of J' turns out to be 
sufficient. 

Determination of and the size dependence of A 
upon L, Al, is essential to the computation of Qa- The 
measurements of the correlation length must be taken 
into consideration for the computation of the scaling vari- 
able x. The correlation lengths are measured employ- 
ing the-standard formula of second moment correlation 
lengthily. For the simplicity of our presentation, here we 
mainly focus on the correlation length. For each J', we 
chose three different inverse temperatures for the com- 
putations of the scaling variable and scaling function, ft 
is important to select the values of the temperatures so 
that they are in the scaling regime. The criterion for 
this, which is known from the studies of well-known pure 
systems such as 2D Ising and three-state Potts models, is 
that the thermodynamic value of the correlation length 
at a temperature in the scaling regime should be suffi- 
ciently large, namely £ > 5. 



IV. RESULT AND ANALYSIS 

The Q c O) are calculated for J' = 0.9,0.25, and 0.1. 
For an illustration, the size dependence of the £ and Xi 
and the results of the computations of x, Q^(x), and 
Q x {x) obtained for J' = 0.25 and (3 — 0.77 and for 
J' = 0.1 and (3 = 1.04 are tabulated in Table (|). It 
is worth noticing that both £l and xl are monotonically 
increasing function of L. (This is the case even in terms 
of the scaling variable s.) The L dependence becomes 
weaker with increasing L, and becomes vanishingly small 
for sufficiently large L. This L independent value within 
the statistical errors is the corresponding thermodynamic 
value. In Table(|) it is shown that one needs L/£l ^ 10 
for the accurate measurements of the thermodynamic val- 
ues. The value of the ratio (L to £l) is larger than wtthe 
corresponding pure system that is approximately 7li3. It 
is worth stressing that the value of the ratio beyond which 
a measurement becomes thermodynamic is temperature 
independent due to Eq. (|l3|) . One thus needs to simulate 
at least with L larger than 10 times of the correlation 
length to make the data be free of finite size effect. 

The plots of Q^(x) for the three values of (3 in the scal- 
ing regime are shown in Fig. (|l|) for J' — 0.25. Clearly 
each data set belonging to different (3 collapses onto a sin- 
gle curve that is the universal FSS function for the value 
of J'. Thus the validity of the FSS is verified for the J'. 



We repeated similar procedure for the other two values of 
J' and observed similar data-collapse. For comparison, 
the FSS function for each J' is shown in our Fig.(||). It 
is observed that the FSS function for J' — 0.9 is indis- 
tinguishable from that of the corresponding pure system 
that was calculated in Ref.( [n]). However, for stronger 
disorder the FSS scaling function clearly depends on J'. 
The ansatz for our scaling function is 



1 + ClX + C 2 X + C 3 X 



CaX 



(18) 



The values of the coefficients ct(i = 1, ... ,4) are calcu- 
lated by fitting our data to the ansatz for each value of 
J'. 

The thermodynamic correlation lengths are measured 
directly (without using any extrapolation method) for 
the data roughly over the range £ < 40. With the 
knowledge of the scaling function available, the thermo- 
dynamic values closer to criticality can now be easily es- 
timated-by the use of the single-step FSS extrapolation 
methodliS. For the each value of J' the thermodynamic 
values of the correlation length are evaluated over the 
ranges 5.7(1) < £ < 204(2) ([5.7(1), 204(2)]), [5.8(1), 
217(3)], and [5.0(1), 203(5)]. For the J' = 0.05 the ther- 
modynamic data were obtained by direct measurement 
under the thermodynamic condition > 12. which 

is over the range 5.52(2) < £, < 47.64(29). Our thermo- 
dynamic data are summerized in Table (O). 

ln£(i) as a function of | lnf| is plotted in Fig.(||). The 
slope of each straight line corresponds to the value of v. It 
is evident that v increases with decreasing J' at least ef- 
fectively. Fixing the critical points at the self-dual points 
in the x 2 fits and assuming a pure power-law type critical 
behavior, we obtain v = 1.01(1), 1.10(2), 1.20(3), 1.34(6) 
for J' = 0.9, 0.25, 0.1, and 0.05 respectively. Assuming 
a scaling function with a nonconfluent correction term, 
e.g., £(£) ~ t~ v (l + at), yields the estimate of the critical 
exponent, e.g., v = 1.08(4) and v = 1.17(5) for J 1 = 0.25 
and J' = 0.1 respectively. Notice that for J' — 0.9 the 
estimated value of v is virtually the same as that in the 
pure system. 

Although our thermodynamic data are extraordinarily 
broad compared to standard Monte Carlo or series ex- 
pansion study can usually make available, it turns out 
that the data fit to a pure power-law as well as to the 
form with the multiplicative logarithmic correction. This 
can be easily seen by simple numerical experiment as 
well. Namely, generate some numbers according to a cer- 
tain power-law with the value of v slightly larger than 
1, and fit them to a multiplicative correction of the form 
Eq. (||) . One can always find acceptable values of C and v 
over quite broad ranges, making it prohibitly difficult to 
get accurate estimate of the logarithmic exponent. With 
the use of the asymptotic form,Eq.(|k can one get much 
shaper estimate of v than with Eq.(5j). The problem is 
that one does not know a priori whether all the data are 
beyond the crossover point. Nevertheless, we find that 
use of the asymptotic form rather than Eq.(|5|) is quite 
useful for our purpose. 
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A useful observation is that [1 + C\ ln(t) |] p (y > 0) is 
less singular than | \nt\ u for any range of t > 0. Namely, 
note that 



(1 + C\ \nt\f ~ | lntT (C + 1/| \nt\Y 



(19) 



and that the second term on the right hand side of 
Eq. ( |l9|) suppresses the singularity for any physical value 
of t. Thus the asymptotic singularity Eq.(||) is al- 
ways more singular than the mixture of the singularities, 
Eq. ([5]) . Accordingly, if a thermodynamic correlation data 
turn out to be more singular than the asymptotic singu- 
larity in an arbitrary portion of the scaling regime, then 
the the prediction of SSL must be invalidated. 

In Fig.(|) we plotted 1 lnt| ' 5 ) for J' = 0.25, 

0.1, and 0.05. For the J' — 0.25, it is observed that the 
value of the ratio decreases monotonically until tempera- 
ture is very close to criticality, but starts to increase with 
further approaching to the criticality. This is surprising 
in view of the SSL's picture, because the figure shows 
none of the data are either in the asymptotic regime or 
in the scaling regime of the pure system. For the J' = 0.1, 
we find that the data less close to the criticality are con- 
sistent with the prediction of SSL, but start to deviate 
from it as t — > 0. For the J 1 = 0.05, we observe that 
the data scale faster than the asymptotic form for all our 
data. We thus lead to the picture that seemingly consis- 
tency with the logarithmic correction for weak disorder 
starts to become invalidated in sufficiently deep scaling 
regime. For the very strong disorder, any portion of data 
are inconsistent with the asymptotic form. 

Binder's cumulant ratio at piticality, denoted by 
is another universal quant ityE3E3. For each J' we mea- 
sured it at the critical point with varying L (TablepTl|) . It 

is observed that given J' is invariant with L within 
the statistical errors, and that it tend to increase uni- 
formly with decreasing J'. The value for J' = 0.9 is 
indistinguishable from the pure case, as is the case for 
the scaling function Q^. 



: random-coupling 



numerical evidence was obtained for t ; 
three-state Potts ferromagnet as wellEiJ. 

As mentioned in the introduction most of previous nu- 
merical studies claiming for the evidence for the SSL 
came from the observation of the consistency of the log- 
arithmic correction. The consistency, however, cannot 
disprove the pure power-law singularity that is essen- 
tially more steeper than the logarithmic correction. A 
very strong claim for the evidence for the SSL was made 
in a recent high temperature expansion study of the same 
modeltij. What the authors of the paper observe is the 
monotonic increase of the effective value of 7 with the 
strength of disorder. On the other hand, they claimed 
that the same data fitted to the prediction of SSL, Eq([|), 
give rise to the value of the logarithmic exponent inde- 
pendent of the disorder. This is unlikely to be mathe- 
matically correct because the effective increase of 7 leads 
to the effective increase of the logarithmic exponent, _as 
can be easily checked by simple numerical experimcntLil. 

The-. s±.uj±jz_. of the FSS behavior of the specific 
heatclo'LJ'LJLJ, instead of its thermodynamic behavior, 
should also be treated with caution. In Ref.( [l0]) it was 
shown that the size dependence of the specific heat is 
different from that of the correlation length or the mag- 
netic susceptibility. Namely, the specific heat increases 
with s = L/£ (t) for smaller s but starts to decreases with 
it for larger s. In view of FSS, Eq.(ll), what this implies 
is the existence of the upper bound value of the scaling 
variable s below which the specific heat increases with L. 
Thus, for any value of L, one is bound to observe mono- 
tonic increase of the specific heat at criticality where £ di- 
verges (in other words, s is smaller than the upper bound 
value for any finite L when £ — > 00), and this has nothing 
to do with the actual divergence of the specific heat at 
criticality. ftom the study of the FSS behavior of the 
pure systemEa, we also note that FSS of the specific heat 
gives rise to less accurate result than that of £ or \- 



V. CONCLUSION AND DISCUSSION 



We have obtained clean numerical evidence that FSS 
holds for quenched random disordered Ising ferromag- 
net. The universal FSS function is found to be depen- 
dent upon the strength of disorder for strongly disordered 
cases. For very weak disorder, however, it appears that 
disorder does not induce new universality class to the res- 
olution of our data. The behavior of Binder's cumulant 
ratio is in complete agreement with the feature obtained 
from the FSS scaling function. We also show that the 
singularity of £ is steeper than the theoretical prediction 
for the data sufficiently close to criticality, which almost 
certainly rules out the validity of SSL in the strong dis- 
ordered system. Our result of varying exponent v com- 
bined with the established fact of the invariance of "//is 
supports the scenario of weak universalityES. The same 



5 



TABLE I. Size dependence of various physical quantities 
and the computed scaling variable and scaling functions at 
/3 = 0.77 with J' = 0.25 (the upper part) and at (3 = 1.04 
with J' = 0.1 (the lower part). The L independent value 
(within the statistical errors) is the corresponding thermody- 
namic value, which are £ = 18.9(2) and £ = 25.4(3) respec- 
tively for 13 = 0.77 and 1.04. 



L 



XL 



Qe 



Qx 



20 
24 
30 
34 
40 
50 
60 
70 
80 
150 
200 
250 



11.51(9) 
12.69(9) 
14.21(12) 
15.02(10) 
15.77(8) 
16.94(10) 
17.42(10) 
17.78(13) 
18.20(18) 
18.71(14) 
18.86(11) 
18.85(17) 



142.0(8) 
178.0(1.2 
231.1(1.7 
262.4(1.7 
301.5(1.8 
356.3(2.9 
390.1(3.3 
408.6(4.0 
428.4(5.6 
456.0(5.0 
456.3(2.5 
456.2(3.8 



0.576(5 
0.529(4 
0.474(4 
0.442(3 
0.394(2 
0.339(2 
0.290(2 
0.254(2 
0.228(2 
0.125(1 
0.094(1 
0.075(1 



0.609(11 
0.671(12 
0.752(14 
0.794(14 
0.834(13 
0.896(15 
0.922(15 
0.941(17 
0.963(20 
0.992(18 
0.998(16 
0.997(21 



20 
30 
50 
80 
120 
160 
200 
240 
300 
360 



12.61(9) 
16.17(14) 
20.41(23) 
23.01(19) 
24.16(17) 
24.55(15) 
24.86(14) 
25.17(14) 
25.35(16) 
25.29(29) 



154.0(8) 
267.1(2.1 
467.7(5.8 
647.1(7.2 
720.1(7.5 
746.5(6.4 
759.0(5.9 
777.7(4.8 
776.2(4.5 
775.3(7.8 



0.630(5 
0.539(5 
0.408(5 
0.288(2 
0.201(1 
0.153(1 
0.124(1 
0.105(1 
0.084(1 
0.070(1 



0.498(7) 
0.639(10 
0.807(15 
0.909(15 
0.955(14 
0.970(14 
0.983(13 
0.995(13 
1.002(14 
1.000(19 



G 



TABLE II. Selection of thermodynamic £ for J'=0.9, 0.25, 0.1 and 0.05. For the J'=0.05 all the results are obtained by 
direct measurements (without the use of the FSS extrapolation method). 



J' 


= 0.9 


P 


0.42 


0.44 


0.45 


0.455 


0.46 


0.462 


0.463 








5.76(5) 


10.7(1) 


18.4(1) 


28.3(2) 


61.9(3) 


116(1) 


204(2) 


J' 


= 0.25 


P 


0.70 


0.74 


0.77 


0.78 


0.79 


0.80 


0.803 








5.80(5) 


9.72(7) 


18.9(1) 


26.7(2) 


43.6(2) 


116(1) 


217(3) 


J' 


= 0.1 


P 


0.87 


0.92 


1.0 


1.04 


1.06 


1.08 


1.093 








5.05(7) 


6.95(8) 


14.20(10) 


25.4(3) 


39.4(4) 


80.5(1.4) 


203(5) 


./' 


= 0.05 


P 


1.00 


1.10 


1.20 


1.25 


1.28 












5.52(2) 


9.15(9) 


18.40(11) 


30.75(18) 


47.64(29) 







7 



TABLE III. Binder's cumulant ratio at the self-dual points, 
for the three values of J'. Note that f7i(i = 0) for each J' does 
not vary with L within the statistical errors, thus showing 
that each self-dual point is indeed the critical point. Also, it 
is clear that Ul (t = 0) increases with decreasing J' , although 
for J' = 0.9 it is hardly distinguishable from the value of the 
pure system. For J' — 0.25 we extended the measurements 
up to L=400, which does not show any sign of crossover. 



L 


J' = 1.00 


J' = 0.90 


J' = 0.25 


J' = 0.10 


20 


1.8324(6) 


1.834(1) 


1.850(3) 


1.862(3) 


40 


1.8321(6) 


1.833(2) 


1.846(3) 


1.858(3) 


60 


1.8317(5) 


1.832(1) 


1.852(3) 


1.854(4) 


80 


1.8318(5) 


1.833(1) 


1.847(3) 


1.862(4) 


100 


1.8316(6) 


1.832(2) 


1.849(3) 


1.858(4) 


200 






1.844(3) 




400 






1.845(2) 
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FIG. 1. Numerical calculation of for J' = 0.25, at three 
arbitrary inverse temperature /3=0.74, 0.77, and 0.79 in the 
scaling regime. Our figure clearly indicates that the scaling 
function does not have explicit temperature dependence. For 
all the values of /3, the smallest value of L used is 20. 

FIG. 2. The Q e (x) for the three values of J'. The depen- 
dence of the scaling function on J' is obvious. 

FIG. 3. ln£ versus |lnt|. The dotted lines represent the 
results of the best x 2 fits assuming pure power-law type sin- 
gularity. The values of the slope , which correspond to the 
values of u, are 1.00, 1.10, 1.20, and 1.34 respectively for 
J' = 0.90, 0.25, 0.10, and 0.05. 

FIG. 4. The ratio $(t) /(t" 1 ! lni| 1/2 ) for J' = 0.25, 0.1, and 
0.05. The data for J' = 0.1 and 0.05 are uniformly shifted so 
that the difference in the data points are clearly visible in a 
single figure. Here increasing value of the ratio as t becomes 
smaller is an evidence that the data over the regime are in- 
consistent with the prediction of SSL. The tendency of the 
increasing ratio with t becoming sufficiently small is observed 
for J' — 0.25 and 0.1, showing that the seemingly consistency 
with the logarithmic correction of SSL for the larger values of 
t becomes invalidated in sufficiently deep scaling regime. For 
the J' — 0.05, all the data scales steeper than the asymptotic 
scaling, showing that SSL cannot be correct for any data in 
the scaling regime for this J'. 
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